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ABSTRACT 

r-j . Foreground removal is a major challenge for detecting the redshifted 21-cm neutral 

Q_j hydrogen (HI) signal from the Epoch of Reionization (EoR). We have used 150 MHz 

GMRT observations to characterize the statistical properties of the foregrounds in four 
different fields of view. The measured multi-frequency angular power spectrum Ct {Ai>) 
is found to have values in the range 10 4 mK 2 to 2 x 10 4 mK 2 across 700 < £ < 2 x 10 4 
j^j , and Av < 2.5 MHz, which is consistent with model predictions where point sources 

are the most dominant foreground component. The measured Ct(Av) does not show 
a smooth Av dependence, which poses a severe difficulty for foreground removal using 
polynomial fitting. 

The observational data was used to assess point source subtraction. Considering 
the brightest source (~ 1 Jy) in each field, we find that the residual artifacts are 
\q . less than 1.5% in the most sensitive field (FIELD I). Considering all the sources in 

the fields, we find that the bulk of the image is free of artifacts, the artifacts being 
localized to the vicinity of the brightest sources. We have used FIELD I, which has a 
rms noise of 1.3 mJy/Beam, to study the properties of the radio source population to 
a limiting flux of 9 mJy. The differential source count is well fitted with a single power 
law of slope —1.6. We find there is no evidence for flattening of the source counts 
towards lower flux densities which suggests that source population is dominated by 
the classical radio-loud Active Galactic Nucleus ( AGN) . 

The diffuse Galactic emission is revealed after the point sources are subtracted 
out from FIELD I . We find C e oc £ -2 ' 34 for 253 < I < 800 which is characteristic 
C$ \ of the Galactic synchrotron radiation measured at higher frequencies and larger an- 

gular scales. We estimate the fluctuations in the Galactic synchrotron emission to be 
y/£(e+l)Ce/2iv ~ 10 K at £ = 800 (9 > 10'). The measured C t is dominated by the 
residual point sources and artifacts at smaller angular scales where Ct ~ 10 3 mK 2 for 
£ > 800. 

Key words: techniques:interferometric-radio continuum:general- (cosmology:)diffuse 
radiation 
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1 INTRODUCTION 

Observations of the bigb-redshift Universe with the 21-cm 
hyperfine line emitted by neutral hydrogen gas (HI) dur- 
ing the Epoch of Reionization (EoR) contains a wealth of 
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cosmological and astrophysical informations (see Furlanetto 
,Oh & Briggs 2006 ; Morales & Wyithe 2010 for recent re- 
views). Analysis of quasar absorption spectra (Becker et al. 
2001; Fan et al. 2002) and the CMBR (Komatsu et al. 2009) 
together indicate that reionization probably started around 
a redshift of 15 and lasted until a redshift of 6. The Gi- 
ant Metrewave Radio Telescope (GMRT ; Swamp et al. 



http://www.gmrt.ncra.tifr.res.in 



2 A. Ghosh, J. Prasad, S. Bharadwaj, S. S. Ali and J. N. Chengalur 



1991) is currently functioning at several bands in the fre- 
quency range of 150 MHz to 1420 MHz and can potentially 
detect the 21 cm signal at high redshifts. Several upcom- 
ing low- frequency telescopes such as LOFAR 2 , MWA 3 and 
PAPER 4 are also targeting a detection of the high redshift 
21-cm signal. LOFAR is already operational and will soon 
start to collect data for the EoR project. 

In this paper we have carried out GMRT observations 
to characterize the statistical properties of the background 
radiation at 150 MHz. These observations cover a frequency 
band of ~ 5 MHz with a ~ 100 kHz resolution, and cover 
the angular scales ~ 30' to ~ 30". The 21 cm radiation is ex- 
pected to have an rms brightness temperature of a few mK 
on angular scales of a few arc minute (Zaldarriaga, Furlan- 
etto & Hernquist 2004). This signal, is however, buried in 
the emission from other astrophysical sources which are col- 
lectively referred to as foregrounds. Foreground models (Ali, 
Bharadwaj & Chengalur 2008) predict that they are dom- 
inated by the extragalactic radio sources and the diffuse 
synchrotron radiation from our own Galaxy which makes 
a smaller contribution. Analytic estimates of the HI 21-cm 
signal indicate that the predicted signal decorrelates rapidly 
with increasing Av and the signal falls off 90% or more 
at Av > 0.5 MHz (Bharadwaj and Sethi 2001; Bharad- 
waj and Ali 2005). The foreground from different astro- 
physical sources are expected to be correlated over a large 
frequency frequency separation. This property of the signal 
holds the promise of allowing us to separate the signal from 
the foregrounds (Ghosh et al. 2011b). In this paper we have 
characterized the foreground properties across a frequency 
separation Av of 2.5 MHz, the signal may be safely assumed 
to have decorrelated well within this frequency range. 

There currently exist several surveys which cover large 
regions of the sky at around 150 MHz like the 3CR sur- 
vey (Edge et al. 1959; Bennett 1962), the 6C survey (Hales, 
Baldwin & Warner 1988; Waldram 1998) and the 7C survey 
(Hales et al. 2007). These surveys have a relatively poor an- 
gular resolution and a limiting flux density of ~ 100 mjy. 
Di Matteo et al. (2002) have used the 6C survey, the 3CR 
survey and the 3 CRR catalogue (Laing , Riley & Longair 
1983) to estimate the extragalactic point source contribution 
at 150 MHz. The GMRT, which is currently the most sensi- 
tive telescope at 150 MHz, offers the possibility to perform 
deeper surveys with a higher angular resolution. A few sin- 
gle pointing 150 MHz surveys (Sirothia et al. 2009; Ishwara- 
Chandra et al. 2010, 2011; Interna et al. 2011) have been 
performed with an rms noise of ~ 1 mJy/Beam at a reso- 
lution of 20" - 30". The TGSS 5 is currently underway to 
survey the sky north of the declination —30° with an rms 
noise of 7 — 9 mJy/Beam at an angular resolution of 20". 

Diffuse synchrotron radiation produced by cosmic ray 
electrons propagating in the magnetic field of our Galaxy 
(Ginzburg & Syrovatskii 1969) is a major foreground compo- 
nent. This component, which is associated with our Galactic 
disk, is particularly strong in the plane of our Galaxy. Our 
target fields were selected well off the Galactic plane with 
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b > 10° where we expect a relatively lower contribution from 
the Galactic synchrotron radiation. 

Radio surveys at 408 MHz (Haslam et al. 1982), 1.42 
GHz (Reich 1982; Reich & Reich 1988) and 2.326 GHz 
(Jonas et al. 1998) have measured the diffuse Galactic syn- 
chrotron radiation at angular scales larger than ~ 1° where 
it is the most dominant contribution, exceeding the point 
sources. Platania et al. (1998) have analyzed these observa- 
tions to show that the Galactic synchrotron emission has a 
steep spectral index of a « 2.8. The angular power spectrum 
Ci has been measured at 2.3 GHz (Giardino et al. 2001) and 
2.4 GHz (Giardino et al. 2002) where they find Ce ~ r 2A at 
angular multipoles £ < 250. WMAP data shows Ce ~ l~ 2 at 
I < 200 (Bennett et al. 2003) which is slightly flatter com- 
pared to the findings at lower frequencies. It is relevant to 
note that all of these results are restricted to angular scales 
greater than 43' (£ < 250), little is known (except Bernardi 
et al. 2009 discussed later) about the structure of the Galac- 
tic synchrotron radiation at the sub-degree angular scales 
probed in this paper. 

Ali, Bharadwaj & Chengalur (2008) have carried out 
GMRT observations to characterize the foregrounds on 
sub-degree angular scales at 150 MHz. They have used 
the measured visibilities to directly determine the multi- 
frequency angular power spectrum Ce(Av) (MAPS; Datta, 
Roy Choudhury & Bharadwaj (2007)) which characterizes 
the statistical properties of the fluctuations in the back- 
ground radiation jointly as a function of the angular multi- 
pole I and frequency separation Av. They find that the mea- 
sured Ce(Av) has a value of around 10 4 mK 2 which is seven 
orders of magnitude larger than the expected HI signal. The 
measured Ct(Av) was also found to exhibit relatively large 
(~ 10 — 20%) fluctuations in Av which pose a severe prob- 
lem for foreground removal. Further, it was attempted to 
remove the foregrounds by subtracting out all the identifi- 
able point sources above 8mJy. The resulting Ci, however, 
only dropped to ~ 3 times the original value due to residual 
artifacts. 

Bernardi et al. (2009) have analyzed 150 MHz WSRT 
observations near the Galactic plane ((I = 137°,?) = +8°) to 
characterize the statistical properties of the diffuse Galac- 
tic emission (after removing the point sources) at angular 
multipoles £ < 3000 (0 > 3.6'). They find that the mea- 
sured total intensity angular power spectrum shows a power 
law behaviour Ce oc £~ 2 ' 2 for £ < 900. They also measured 
the polarization angular power spectrum for which they find 
Ce oc l- 1 - 65 for £ < 2700. Their observations indicate fluc- 
tuations of the order of ~ 5.7 K and 3.3 K on 5' angular 
scales for the total intensity and the polarized diffuse Galac- 
tic emission respectively. Pen et al. (2009) have carried out 
150 MHz GMRT observati ons at a h igh Galactic latitude to 
place an upper limit of y/PCe/2-K < 3K on the polarized 
foregrounds at I < 1000. 

There has been a considerable amount of work towards 
simulating the foregrounds (Jelic et al. 2008; Sun et al. 2008; 
Bowman, Morales & Hewitt 2009; Gleser et al. 2008; Harker 
et al. 2009; Liu, Tegmark & Zaldarriaga 2009; Liu et al. 
2009; Pctrovic & Oh 2011) in order to develop algorithms 
to subtract the foregrounds and detect the redshifted 21- 
cm signal. These simulations require guidance from obser- 
vational data, and it is crucial to accurately characterize the 
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foregrounds in order to have realistic simulations for future 
observations. 



One of the challenges for the detection of this cosmo- 
logical 21-cm signal is the accuracy of the foreground source 
removal. Our current paper aims at characterizing the fore- 
ground at arcminute angular scales which will be essential 
for extracting the 21-cm signal from the data. In this paper 
we study the radio source population at 150 MHz, which 
is still less explored, to determine the differential source 
count which plays a very important role in predicting the 
angular power spectrum of point sources (Ali, Bharadwaj 
& Chengalur 2008; Ghosh et al. 2011a). We find that the 
point sources are the most dominating foreground compo- 
nent at the angular scales of our analysis. The accuracy of 
point source removal is one of the principal challenges for 
the detection of the 21-cm signal. We investigate in detail 
the extragalactic point source contamination in our observed 
fields and study how accurately the bright point sources can 
be removed from our images. We also estimate the level of 
residual contamination in our images. It is expected that 
the diffuse Galactic emission will be revealed if the point 
sources are individually modelled and removed with high 
level of precession. To our knowledge, the present work is 
only the second time that the fluctuations in the Galactic 
diffuse emission have been detected at around 10 arcminute 
angular scales at a frequency of 150 MHz . We also comple- 
ment this with a measurement of the angular power spec- 
trum of the diffuse Galactic foreground emission in the £ 
range 200 to 800. 



In this paper we have analyzed GMRT observations 
to characterize the foregrounds in four different fields 
at 150 MHz which corresponds to the HI signal from 
z = 8.3. We note that unless otherwise stated we use 
the cosmological parameters (f2 m o, Oao, Qbh 2 , h, n s , os) = 
(0.3, 0.7, 0.024, 0.7, 1.0, 1.0) in our analysis. We now present 
a brief outline of the paper. Section 2 describes the observa- 
tions and data analysis. We have used the multi-frequency 
angular power spectrum Ce(Au) to quantify the statisti- 
cal properties of the background radiation in the observed 
fields. Section 3 presents the technique that we have used 
to estimate C'i(Av) directly from the measured visibilities. 
We show the the measured Ce{Av), and discuss its prop- 
erties in Section 4. In this section we have also compared 
the measured Ci(Av) with foreground model predictions. In 
Section 5 we consider the extragalactic point sources which 
forms the most dominant foreground component in our ob- 
servations. In order to assess how well we can subtract out 
the point sources, we consider point source subtraction for 
the brightest source in each field, as well as subtracting out 
all the identifiable sources. We use the most sensitive field 
(FIELD I) in our observation to study the nature of the 
radio source population at 150 MHz, and we also provide a 
source catalogue for this field in the online version of this 
paper (Appendix A) . In Section 6 we have used the residual 
data of FIELD I, after point source removal, to study the 
diffuse Galactic synchrotron radiation. We have summarized 
the results of the entire paper in Section 7. 



2 GMRT OBSERVATIONS AND DATA 
ANALYSIS 

The GMRT has a hybrid configuration (Swarup et al. 1991), 
each of diameter 45 m, where 14 of the 30 antennas arc 
randomly distributed in a Central Square which is approxi- 
mately 1.1km x 1.1km in extent. The rest of the antennas 
lie along three ~ 14 km long arms in an approximately 'Y' 
configuration. The shortest antenna separation (baseline) is 
around 60 m including projection effects while the largest 
separation is around 26 km. The hybrid configuration of the 
GMRT gives reasonably good sensitivity to probe both com- 
pact and extended sources. Table 1 summarizes our obser- 
vations including the calibrators used for each field. For all 
observations, visibilities were recorded for two circular po- 
larisations (RR and LL) with 128 frequency channels and 
16 s integration time. 

We have observed FIELD I in GTAC (GMRT Time 
Allocation Committee) cycle 15 in January 2008, whereas 
FIELD II and III were observed in cycle 17 during Febru- 
ary 2010. Our target fields were selected at high Galactic 
latitudes (b > 10°) which were up at night time during the 
GTAC cycle 15 and 17, and which contain relatively few 
bright sources (> 0.3 Jy) in the 1400 MHz NVSS survey. 
We have selected such fields where the 408 MHz Haslam 
map sky temperature is relatively low (in the range of 30 K 
to 40 K) and there is no significant structure visible at the 
angular resolution (~ 0.85°) of Haslam map. FIELD IV had 
been observed in cycle 8 (June 2005) and the 150 MHz fore- 
grounds in this field have been analyzed and published ear- 
lier (Ali, Bharadwaj & Chengalur 2008). We have used the 
calibrated data set of Ali, Bharadwaj & Chengalur (2008), 
the foreground characterisation has, however, been redone 
using an improved technique (Ghosh et al. 2011b). 

The observations were largely carried out at night to 
minimise the Radio Frequency Interference (RFI) from man- 
made sources. Further, the ionosphere is considerably more 
stable at night, and the scintillation due to ionospheric os- 
cillations is expected to go down. The duration of our ob- 
servations (~ 4 - 10 hrs.) was chosen so as to achieve a 
reasonably high signal to noise ratio, and also to have a rea- 
sonably large number of visibilities adequate for a statistical 
analysis of the foregrounds. 

We have followed the standard procedure of interfero- 
metric observations where we observe a flux calibrator of 
known strength and a phase calibrator whose structure is 
known but whose strength varies from epoch to epoch. The 
flux calibrator, which is observed at the start of the observa- 
tion, sets the amplitude scale of the gains and this is used to 
set the amplitude of the phase calibrator. Phase calibrators 
were chosen near the target fields, and these were observed 
every half an hour so as to correct for temporal variations in 
the system gain. We have used the observations of the phase 
calibrators to determine the complex gains of the antennas, 
and these gains were interpolated to the observations of the 
target fields. 



The data for FIELDS I, II and III were reduced us- 
ing FLAGCAL (Prasad & Chengalur (2012) which is a 
flagging and calibration software for radio interferometric 
data. FIELD IV, however, was not analyzed with FLAG- 
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Table 1. Observation summary 





FIELD I 


FIELD II 


FIELD III 


FIELD IV 


Date 


2008 Jan 08 


2010 Feb 07 


2010 Feb 08 


2005 June 15 


Working antennas 


28 


28 


26 


23 


Central Frequency 


153 MHz 


148 MHz 


148 MHz 


153 MHz 


Channel width 


62.5 kHz 


125.0 kHz 


125.0 kHz 


62.5 kHz 


Bandwidth 


3.75 MHz 


8.75 MHz 


10.0 MHz 


4.37 MHz 


Total observation time 


11 hrs 


6 hrs 


6 hrs 


13 hrs 


Flux calibrator 

Observation time 

Flux density (Pcrley- Taylor 99) 


3C147 
1.2 hrs 
~ 72.5 Jy 


3C147 
1.75 hrs 
~ 72.6 Jy 


3C286 
0.5 hrs 
~ 33.41 Jy 


3C48 
2 hrs 
~ 62.5 Jy 


Phase calibrator 
Observation time 
Flux density 


3C147 


3C147 


1459+716 

1.6 hrs 
~ 38.7 Jy 


3C48 


Target field (a, <5) 2 ooo 

Galactic coordinates (/, b) 

Sky Temp. (Haslam Map, 408 MHz) 

Observation time 


(05 h 30 m 00 s , 
+60°00'00") 
151.80°, 13.89° 
~ 40 K 
9.8 hrs 


(06 h 00 m 00 s , 
+62° 12' 58") 
151.49°, 18.12° 
~ 30 K 
4.25 hrs 


(12 ,I 36 m 49 s , 
+62° 12' 58" ) 
125.89°, 54.83° 

~ 20 K 

3.9 hrs 


(l h 36 m 48 s , 
+41°24'23") 
132.00°, 20.67° 
~ 30 K 
11 hrs 



Table 2. Wide field Imaging summary 





FIELD I 


FIELD II 


FIELD III 


FIELD IV 


Image size 


5700 x 5700 


5500 x 5500 


5500 x 5500 


5500 x 5500 


Pixel size 


3.22" x 3.22" 


3.32" x 3.32" 


3.27" x 3.27" 


3.40" x 3.40" 


Number of facets 


121 with 2 overlap 


109 with 2 overlap 


109 with 2 overlap 


109 with 2 overlap 


Off-source Noise 


1.3 mJy/Beam 


2.5 mJy/Beam 


4.5mJy/Beam 


1.6mJy/Beam 


Peak/Noise 


700 


620 


422 


550 


Flux density (max., min.) 


(905 mJy/Beam, 
— 14 mJy/Beam) 


(1.55 Jy/Bcam, 
-28 mJy/Beam) 


(1.9 Jy /Beam, 
—45 m Jy/Bcam) 


(900 mjy /Beam, 
— 47mJy/Beam) 


Synthesized beam 


21" x 18", PA=61° 


30" x 22", PA=-52° 


33" x 20", PA=53° 


24" x 18", PA=70° 



CAL, it was flagged manually and then calibrated using the 
standard task within Astronomical Image Processing Soft- 
ware (AIPS). The details of the observation, data reduction 
and relevant statistics for FIELD IV are discussed in Ali, 
Bharadwaj & Chengalur (2008). 

RFI is a major factor limiting the sensitivity of at low 
frequencies. RFIs effectively increases the system noise and 
corrupts the calibration solutions. It also restricts the avail- 
able frequency bandwidth. The effect is particularly strong 
at frequencies below 0.5 GHz. FLAGCAL identifies and re- 
moves bad visibilities by requiring that good visibilities be 
continuous in time and frequency, computes calibration so- 
lutions using known flux and phase calibrators and interpo- 
lates these onto the target fields (see Prasad & Chengalur 
(2012) for more details ). The flagged and calibrated visi- 



bility data were used to make continuum images using the 
standard tasks in the AIPS. 

The large field of view (hereafter FoV; #fwhm = 3.8°) 
of the GMRT at 150 MHz leads to considerable error if the 
non-coplanar nature of the GMRT antenna distribution is 
not taken into account. We use the three dimensional (3D) 
imaging feature ( Perley 1999) in the AIPS task 'IMAGR' 
in which the entire field of view is divided into multiple 
sub-fields (facets) , each of which is imaged separately. Table 
2 contains a summary of the imaging details with all the 
relevant parameters which are mostly self-explanatory. 

The presence of a large number of bright sources in the 
fields that we have observed allows us to carry out self cali- 
bration to improve the complex gains. This reduces the er- 
rors from temporal variations in the system gain, and spatial 
and temporal variations in the ionospheric properties. For all 
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the target fields, the data went through 3 rounds of phase 
self calibration followed by one round of self calibration for 
both amplitude and phase. The phase variations often oc- 
cur on timescales of as few minutes or less. The time inter- 
val for gain correction was chosen as 5, 3, 2 and 5 minutes 
for the successive self calibration loops. In successive self- 
calibration loops, the gain solution converged rapidly with 
less than 0.1% bad solutions in each loop. The final gain 
table was applied to all the 128 frequency channels. 

The final continuum images for the four fields are shown 
in Figure 1. The images all have a FoV of 4.0° x4.0° . To avoid 
bandwidth smearing in the continuum image, we have col- 
lapsed adjacent channels within < 0.7 MHz and separately 
combined the respective images. 

The subsequent analysis was done using the final; cali- 
brated visibilities of the original 128 channel data. The final 
data contains visibilities for 2 circular polarisations. The vis- 
ibilities from the two circular polarisations were combined 
(V = [Vrr + Vll]/2) for the subsequent analysis. 



3 ESTIMATION OF ANGULAR POWER 
SPECTRUM FROM VISIBILITY 
CORRELATIONS 

In this paper we have used the multi-frequency angular 
power spectrum Ci(Av) (Datta, Roy Choudhury & Bharad- 
waj 2007) to quantify the statistical properties of the sky 
signal. This jointly characterizes the dependence on angu- 
lar scale i~ x and frequency separation Av. In this section 
we briefly present how we have estimated Ci(Av) from the 
measured visibilities. 

For a frequency v, the angular dependence of the bright- 
ness temperature distribution on the sky T(n, v) may be 
expanded in spherical harmonics as 

ae m (v) Yt m (n) . (1) 
The mult i- frequency angular power spectrum is defined 

as 

C'i {Av) = C'i (v, v + Av) = {a im (v) a* lm {y + Av)) . (2) 

where I refers to the angular multipoles on the sky and 
Av is the frequency separation. 

It is possible to use the correlation between pairs of vis- 
ibilities V(U, v) and V(U + AU, v + Av) to estimate C t (Av) 
where I — 2nU. 

V2(\J,v;U+AU,v+Av) = {V(U,i/)V*(U+AU,i/+Ai/))(3) 

The correlation of a visibility with itself is excluded to 
avoid a positive noise bias in the estimator. Ali, Bharadwaj 
& Chengalur (2008) as well as Dutta et al. (2009) contain 
detailed discussions of the estimator. 

For the sky signal which is dominated by foregrounds 
the measured Ce(Av), for a fixed £, is expected to vary 
smoothly with Av and remain nearly constant over the ob- 
servational bandwidth. In a recent 610 MHz observations 
(Ghosh et al. 2011a) we find instead that in addition to a 
component that exhibits a smooth Av dependence, the mea- 
sured Ci(Av) also has a component that oscillates as a func- 
tion of Av (around 1-4%). We note that similar oscillations, 
with considerably larger amplitudes, have been reported in 



earlier GMRT observations at 153 MHz (Ali, Bharadwaj & 
Chengalur 2008). In our recent work (Ghosh et al. 2011b) 
we have noted that these oscillations could arise from bright 
continuum sources located at large angular separations from 
the phase center. The problem can be mitigated by taper- 
ing the array's sky response with a frequency independent 
window function W(6) that falls off before the first null of 
the primary beam (PB). RFI sources, which are mostly lo- 
cated on the ground, are picked up through the sidelobes of 
the primary beam. Tapering the array's sky response is also 
expected to mitigate the RFI contribution. 

Close to the phase centre, the 150 MHz GMRT PB is 
reasonably well modeled by a Gaussian 
where the parameter 0o is related to the FWHM of the PB 
as 9o « 0.6 x #fwhm, and 6o = 2.3° (#fwhm = 3.8°). We ta- 
per the array's sky response by convolving the observed vis- 
ibilities with a suitably chosen function W(U). The sky re- 
sponse of the convolved visibilities V c — W"®V is modulated 
by the window function W(6) which is the Fourier transform 
of W(XJ). We have used a window function W(6) = e~ e /B ™ 
with W < 9q to taper the sky response so that it falls off 
well before the first null. We parametrise 9 W as 6 W = f6>o 
with f < 1 where do here refers to the value at the fixed 
frequency 150 MHz. 

We have used the convolved visibilities to determine the 
two-visibility correlation which is defined as 

V2(U i ,Ai/) = Vc(U i ,i/)V;(U i ,i/ + Ai/). (4) 

Here U; refers to different points (labelled by i) on 
a grid in the uv plane. We have used a grid of spacing 
AU g = \/ln 2/(ti6 w ) — O.2650™ 1 which corresponds to half 
of the FWHM of W(U), and we have estimated V 2 (Ui, Av) 
at every grid point using all the baselines within a disk of 
radius 2 AU g centered on that grid point. Using eq. (4), how- 
ever, introduces an undesirable positive noise bias in the es- 
timated V2(Uj, Av) (e.g. Begum et al. 2006). It is possible 
to avoid the noise bias if we use the estimator 

V 2 (Vi,Av) = K- 1 x [W(Ui -U a )H>*(U, - U 6 ) 

x V(U a ,z,)V*(U 6 ,^ + A^)] (5) 

where K = £ W(Vi-V a )W*(Vi-V b ) is a normal- 
ization constant and U a , U;, refer to the different baselines 
in the observational data. The noise bias is avoided by drop- 
ping the self-correlations (i.e.. the terms with a — b). We fi- 
nally determine Ci(Av) using (Ali, Bharadwaj & Chengalur 
2008) 

V 2 (U,Av) = ^^y Ce (Av)Q(Av). (6) 

where Q(Av) is a slowly varying function of Av which 
accounts for the fact that we have treated 9 W and {8I v /dT) 
as constants in our analysis. We have used Q(Av) = 1 which 
introduces an extra Av dependence in the estimated Ct(Av) 
. This, we assume, will be a small effect and can be accounted 
for during foreground removal. The gridded data has been 
binned assuming that the the signal is statistically isotropic 
in U. We finally invert equation 6 to determine the angular 
power spectrum 
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Figure 1. The figure shows the continuum images of bandwidth 3.75 MHz, 8.75 MHz, f 0.0 MHz and 4.37 MHz of FIELD I, FIELD II, 
FIELD III and FIELD IV respectively. The off source rms noise is around 1.3mJy/Beam, 2.5mJy/Bcam, 4.5mJy/Bcam and 1.6mJy/Beam 
for FIELD I, FIELD II, FIELD III and FIELD IV respectively. 



where, 6 W is in radians, v is the central observing fre- 
quency in MHz and the Ce(Av) is in mK 2 . 

The measured l/2(U,Ai/) will, in general, have a real 
and imaginary parts (eg. Figure 2). The expectation value of 
V2(U, Av) (eq. 6) is predicted to be real, and the imaginary 
part is predicted to be zero . We use the real part of the 
measured V2(U,Av) to estimate Ct(Av) using eq. (7). A 
small imaginary part, however, arises due to the noise in the 
individual visibilities. The fact that we are correlating two 
slightly different visibilities also makes a small contribution 
to the imaginary part. Figure 2 shows the real and imaginary 
parts of the measured V2(U, Av = 0) as a function of U for 
all the four fields that we have analyzed. As expected, the 



imaginary part is much smaller than the real part. We use 
the requirement that the imaginary part of V2(U, Av) should 
be small compared to the real part as a consistency check of 
our analysis. It also establishes that we are truly measuring 
a sky signal, and our results are not dominated by noise or 
other local effects in the individual visibilities. 

We next investigate whether tapering the sky response 
of the PB actually mitigates the oscillations reported in 
earlier work at 150 MHz (Ali, Bharadwaj & Chengalur 
2008, FIELD IV of this paper). For this purpose, we in- 
troduce the dimensionless decorrelation function Ke(Av) = 
Ct(Av)/Ce(0) which has the maximum value ne(Av) = 1 
at Av = 0, and is in the range | m (Av) |< 1 for other val- 
ues of Av. We use Kt(Av) to quantify the Av dependence 
of Ce(Av) at a fixed value of i. Figure 3 shows ne(Av) for 



Foreground for EoR : 150 MHz GMRT observation 7 



10 
1 

0.1 
0.01 t 
10 t 

1 

0.1 t- 

0.01 



L FIELD I 



•H r- 



\ / 
V 



v I ✓ 



H h 



FIELD IV 



\ / \ 
\ / \ 



11 



10 

1 

0.1 
10 

1 

0.1 



I — h 



i — 

FIELD I 



FIELD 



0.2 



0.5 



0.2 



U kA 



0.5 

U kA 



Figure 2. This shows the real (upper curve) and imaginary (lower curve) parts of the observed visibility correlation V2{U, 0) as a function 
of baselines U. The la error-bars shown on the plot includes contribution from both cosmic variance and system noise. 



FIELD IV with and without the tapering. We have con- 
sidered f = 0.6 and 0.8 which respectively correspond to a 
tapered sky response with FWHM 2.3° and 3.04° as com- 
pared to the untapered PB which has a FWHM of 3.8° at 
150 MHz. The oscillatory pattern is distinctly visible when 
the tapering is not applied. For most values of £, the oscilla- 
tions are considerably reduced and are nearly absent when 
tapering is applied. We do not, however, notice any particu- 
lar qualitative trend with varying f . The tapering, which has 
been implemented through a convolution, is expected to be 
most effective in a situation where the uv plane is densely 
sampled by the baseline distribution. Our results are limited 
by the patchy uv coverage of the observational data that is 
being analyzed here. This also possibly explains why some 
small oscillations persist even after tapering is applied. 

The results for the other fields are very similar to 
FIELD IV, and we have not shown these here. We have 
used a tapered PB with / = 0.8 for all the fields in the 
entire subsequent analysis. 



4 THE MEASURED MAPS C e (Av) 

Figure 4 shows the measured Ci(Av) and 1 — a error es- 
timates for two of our observed fields that we have ana- 
lyzed. We note that the 1 — a error estimate has two dif- 
ferent contributions, the cosmic variance and system noise, 
added in quadrature. In our observations the error is mainly 
dominated by the cosmic variance due to the limited num- 
ber of independent estimates. The system noise makes a 
smaller contribution. We have determined Ct(Av) in the 
range < Av < 2.5 MHz and 700 < I < 2 x 10 4 for FIELDS 
I and IV. FIELDS II and III have been observed for shorter 
time period (Table 1), and the measured Ct{Av) becomes 
relatively noisy at the large baselines which are sparsely 
sampled in our observational data. The £ range has been 
restricted to 700 < t < 10 4 for FIELDS II and III. This 
corresponds to the angular scales 64 to 15 . 

A visual inspection of Figure (4) shows certain regions 
(large i and A^) where we do not have estimates of Ci(Av), 
and these regions have been blanked in the figure. This arises 
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Figure 3. The measured Ki(Av) as a function of Av for the 
different £ values as shown in each panel, for FIELD IV. The thick 
solid, thin solid and dotted curves show results for no tapering, 
and tapering with f = 0.8 and 0.6 respectively. 



due to a combination of two factors. To start with, the large 
baselines are sparsely sampled in the observational data. 
Subsequent flagging of the bad channels further depletes the 
data. These two factors together lead to a situation where 
there is practically no data at large Av separations in the 
large baselines. It is reassuring to note that this effect is very 
similar in all the four FIELDS, of which one was manually 
flagged and the rest underwent automated flagging. This es- 
tablishes that this effect is not an artifact introduced by the 
flagging technique. 

Continuing with the visual inspection of Figure 4, we do 
not find any significant pattern in the variation of Ci(Av) 
with either £ or Av. However, the values seem to change 
more along £ as compared to Av. For each of the four fields, 
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the values of Ce(Au) are nearly all within ±5a of the mean 
value of Ci{Av) which is 5146,9117, 14,000 and 4939 mK 2 
for FIELDS I, II, III and IV respectively. The corresponding 
values of a are shown in the right side of each panel of Figure 
4. The measured Ce(Au) is dominated by the contribution 
from extragalactic point sources, and for each field the mean 
value of Ce(Au) is determined by the flux density of the 



brightest source (S c ) in that field. The variation in the mean 
value of Ce(Au) across the different fields is primarily due 
to the variation in maximum source flux density S c which 
has values S c = 0.905 Jy, 1.55 Jy, 1.9 Jy, 0.9 Jy in FIELDS I, 
II, III and IV respectively (Table 1). 

We next fix Av and study the I dependence of Ct = 
Ci(Av — 0) (Figure 5). The variation in the values of C'e 
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Figure 6. The thick solid line shows the measured Ce as a function of I. The thin solid line with 1 — a error bars shows the total 
foreground contribution based on the model prediction of Ali, Bharadwaj & Chcngalur 2008. The dashed line shows the contribution 
from the Poisson fluctuation of the point sources which is the most dominant contribution at large I. 



do not seem to exhibit any particular pattern, and Ce ap- 
pears to fluctuate randomly around the mean values quoted 
earlier. These fluctuations, as noted earlier, are nearly all 
within ±5cr of the mean values (Figure 5). We see that the 
I dependence of the measured Ce is consistent with random 
fluctuations arising from the statistical uncertainties in the 
measured Ce- Note that the 5 — a error bars shown in figure 
5 reflect the estimated uncertainties arising from both, the 
cosmic variance and the system noise. 

We model the measured Ct using the foreground model 
presented in Ali, Bharadwaj & Chengalur (2008). The model 
predictions, shown in Figure 6, are dominated by the contri- 
bution from extragalactic point sources. At small £, or large 
angular scales, the signal is dominated by the angular clus- 
tering of the extragalactic point sources, whereas the Pois- 
son fluctuation due to the discrete nature of these sources 
is the dominant contribution at large £ which corresponds 
to small angular scales. This transition between these two 
occurs somewhere in the £ range 5,000 to 10 4 , and it is dif- 
ferent in the four fields. The transition shifts to smaller £ 
values in the fields which have a larger value of S c . The 
uncertainty in the model predictions is dominated by the 
Poisson fluctuation of the discrete point sources, and this is 
nearly constant across the entire £ range. The measured Ce 
is consistent with the predictions of the foreground model. 
Nearly all the measured Ce values lie within the la error 
bars of the model predictions. Other astrophysical sources 
like the diffuse synchrotron radiation from our own Galaxy, 
the free-free emissions from our Galaxy and external galax- 
ies (Shaver et al. 1999) make much smaller contributions, 
though each of these is individually larger than the HI sig- 
nal. 

Both simulations (Jelic et al. 2008) and analytical pre- 
dictions (Zaldarriaga, Furlanetto & Hernquist 2004; Bharad- 
waj and Ali 2005; Cooray & Furlanetto 2005; Santos et al. 
2005) suggest that at low £ (£ ~ 10 3 ) the EoR HI signal 
is approximately Cf 1 ~ 10~ 3 — 10~ 4 mK 2 , while at the 
larger £ values (£ ~ 10 4 ) Cf 1 drops to 10~ 5 - 10~ 6 mK 2 . 
We find that the measured Ce, arising from foregrounds, 
has values around 10 4 mK 2 at £ ~ 1000, and it drops by 



50% at the smallest angular scales probed by our observa- 
tions (£ ~ 10 4 ). This suggests that the expected HI signal 
Cf 1 is (~ 10 7 ) times smaller than the measured Ce arising 
from foregrounds. 

We next shift our attention to the Av dependence of 
the measured Ce(Av), holding £ fixed. We have restricted 
our analysis to large angular scales (£ < 2, 000) where the 
HI signal is predicted to be relatively stronger, and we have 
a higher chance of an initial detection. We find that for 
nearly all the £ values for FIELD I and IV the variation in 
Ce(Av) with Av is roughly between 4 x 10 3 mK 2 to 10 4 mK 2 , 
whereas for FIELD II and III Ce(Av) roughly varies between 
6 x 10 3 mK 2 and 3.5 x 10 4 mK 2 across the 2.5 MHz Av range 
that we have analyzed. The fractional variation in Ce(Av) 
ranges from ~ 20 to 40 percent for all the fields. Our fore- 
ground model prediction shows that the measured 150 MHz 
sky signal is dominated by extragalactic point sources. These 
are believed to have slowly varying, power law v a frequency 
spectra. We expect the resultant Ce(Av) to have a smooth 
Av dependence, and essentially show very little variation 
over the small Av range (2.5 MHz) that we have consid- 
ered here (Av/v <C 1). For a fixed £, the cosmic variance is 
expected to introduce the same error (independent of Av) 
across the entire band. As a consequence we do not show 
the cosmic variance for the Av dependence of the measured 
Ce(Av). 

The measured Ce(Av) is, by and large, found to de- 
crease with increasing Av if we hold £ fixed. However, con- 
trary to our expectations, the Av dependence is not com- 
pletely smooth. There are a few £ values, particularly in 
FIELDS II and III, where the Av dependence appears to 
be rather smooth. We however have abrupt variations and 
oscillations in the Av dependence for most of the remain- 
ing data. FIELDS I and IV, which have a channel width 
of 62.5 KHz, have a relatively high frequency resolution as 
compared to FIELDS II and III which have a channel width 
of 125 KHz. For FIELDS I and IV, an oscillatory pattern is 
clearly visible in Ce(Av) at nearly all the £ values. These os- 
cillations, however, are not noticeable in FIELDS II and III 
which have a lower frequency resolution. It is possible that 
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the oscillations are actually also present in FIELDS II and 
III, but we are missing them due to the lower frequency res- 
olution in these two fields. To test this we have also analyzed 
FIELDS I and IV at a lower frequency resolution by collaps- 
ing the original data. We find that the oscillations seen in 
FIELDS I and IV are somewhat reduced, but they can still 
be made out at the lower frequency resolution of 125 KHz. 
We find that though the oscillations in all the fields have re- 
duced considerably after tapering the sky response (Section 
3), a residual oscillatory patterns still persist. The oscilla- 
tions are most pronounced at the lowest i values where the 
period of oscillation varies between ~ 1 MHz — 2 MHz (Fig- 
ures 7) . The period and amplitude of these oscillations both 
decrease with increasing £. The exact cause of this resid- 
ual oscillatory pattern is, at the moment, not known. We 
recollect that the sky tapering is implemented through a 
convolution whose efficiency depends on the uv coverage of 
the baselines. It is possible that the residual oscillations are 
a consequence of the finite and sparse uv coverage of our 
data. 

The oscillations and abrupt changes in the Av depen- 
dence of Ci(Av) posse a severe impediment for foreground 
removal. Foreground removal relies on the key assumption 
that the foregrounds have a Av dependence which is dis- 
tinctly different from that of the HI signal. The foreground 
contributions are expected to vary slowly with increasing 
Av, whereas the HI signal is expected to decorrelate rapidly 
well within Av < 0.5 MHz (Bharadwaj and Sethi 2001; 
Bharadwaj and Ali 2005; Ghosh et al. 2011a). It is then 
possible to remove the foregrounds by modelling and sub- 
tracting out any component that varies slowly with increas- 
ing Av. For example, it is possible to use low order polyno- 
mials to model the Av dependence of the measured Ce(Av) 
at fixed valued of £, and use these to subtract out the fore- 
ground contribution (Ghosh et al. 2011a,b). It is quite evi- 
dent that the oscillations and abrupt changes in the Av de- 
pendence of the measured Ct(Av) posse a severe challenge 
for this technique, or any other technique that is based on 
the assumption that the foregrounds vary smoothly with 
frequency. It is relevant to note that Jelic et al. (2008); 
Harker et al. (2009) have shown that the foreground subtrac- 
tion using polynomial fitting can easily cause over-fitting, in 
which we fit away some of the cosmological signal, or under- 
fitting, in which the residuals of the foreground emission can 
overwhelm the cosmological signal. Therefore, one needs to 
use non-parametric methods that allow data to determinate 
their shape without selecting any a priori functional form of 
the foregrounds (Harker et al. 2009; Chapman et al. 2012). 

We would like to point out that a possible line of ap- 
proach in foreground removal is to represent the sky signal 
as an image cube where in addition to the two angular coor- 
dinates on the sky we have the frequency as the third dimen- 
sion (Jelic et al. 2008; Harker et al. 2009; Chapman et al. 
2012; Zaroubi et al. 2012). For each angular position, polyno- 
mial fitting is used to subtract out the component of the sky 
signal that varies slowly with frequency. The residual sky sig- 
nal is expected to contain only the HI signal and noise (Jelic 
et al. 2008; Bowman, Morales & Hewitt 2009; Liu, Tegmark 
& Zaldarriaga 2009; Harker et al. 2009, 2010). Liu, Tegmark 
& Zaldarriaga (2009) showed that this method of foreground 
removal has problems which could be particularly severe at 
large baselines if the uv sampling is sparse. We would also 
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Figure 7. For FIELD I, this shows the measured C'i(Av) as a 
function of Av for the fixed values of I shown in the panels. The 
error-bars show 10<r system noise. 



like to point out that in a typical GMRT observation the 
uv plane is not completely sampled and due to this limited 
baseline coverage there exist correlated noise in the image 
plane which is a problem for estimating the power spectrum 
in the image plane (for details Dutta (2011), Appendix A). 
Although, the upcoming low frequency radio telescopes such 
as the LOFAR, MWA etc. will have a much dense uv cover- 
age compared to GMRT, where in the core region most of the 
baselines relevant for the EoR experiments will be sampled 
and the foreground removal technique in image-frequency 
space can also be applied (Jelic et al. 2008; Harker et al. 
2009; Chapman et al. 2012; Zaroubi et al. 2012). 



5 POINT SOURCES 

The discrete sources seen in Figure 1 dominate the 150 MHz 
radio sky at the angular scales probed in our observations. 
These sources are typically smaller than the angular resolu- 
tion of our observations. It is now well accepted that these 
are extragalactic, being mainly associated with active galac- 
tic nuclei (AGN). It is of considerable interest to study the 
properties of these sources (Dunlop & Peacock 1990; Jack- 
son & Wall 1999; Peacock 1999; Seymour et al. 2004; Garn 
et al. 2007). 

There currently exist several radio surveys at 150 MHz 
like the 3CR survey (Bennett 1962), the 3CRR catalogue 
(Laing , Riley & Longair 1983), the 6C (Hales, Baldwin & 
Warner 1988; Waldram 1998) and the 7C survey (Hales et 
al. 2007), which cover large regions of the sky. The 6C and 
7C surveys have angular resolutions of ~ 4.2 and ~ 70 
respectively, and the limiting source flux density of these 
surveys is ~ 100 mjy. 

Point sources are the most dominant foreground com- 
ponent and it creates a major problem for detecting the 
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Figure 8. For FIELD II, this shows the measured Ci(Au) as a 
function of Av for the fixed values of I shown in the panels. The 
error-bars show 10cr system noise. 



21-cm signal from the redshifted HI emission in the I range 
10 3 - 10 4 probed here. Di Matteo et al. (2002) have used 
the results of the previous surveys to estimate the point 
source contribution to the foregrounds at 150 MHz. How- 
ever, these have a limiting source flux density of ~ 100 mjy, 
and the extrapolation to fainter sources is rather uncertain. 
The GMRT has an angular resolution of ~ 20" at 150 MHz. 
We are able to achieve a rms noise of around 1.3mJy/Beam 
in the observations reported here. The GMRT is currently 
the only instrument capable of achieving this level of sen- 
sitivity in terms of angular resolution and rms noise. We 
note that there are relatively few GMRT results regarding 
the radio source population at 150 MHz (Ishwara-Chandra 
et al. 2011; Interna et al. 2011) with sensitivity compara- 
ble to that reached in our observation. Recent simulations 
(Bowman, Morales & Hewitt 2009; Liu et al. 2009) also 
indicate that point sources should be subtracted down to a 
~ 10 — 100 mjy threshold in order to detect the EoR signal. 
In this section we use our 150 MHz observations to explore 
the population of radio sources in our observed fields down 
to the detection limit of ~ 9mJy. 

In order to detect the HI signal, it is very important to 
correctly identify the point sources and subtract these out at 
a high level of precision (Bernardi et al. 2011; Pindor et al. 
2011). It is quite evident that a template of detected point 
sources above a given threshold level is an essential part of 
foreground removal where we have noticed that the Poisson 
and clustering component of the point sources are the most 
dominating foreground component at our angular scales of 
analysis. We first focus on the brightest source in each of the 
four fields that we have analyzed. In each field, the brightest 
source (Table 2) alone contributes around 10% of the total 
measured Ct. We consider the brightest source as a test case 
to investigate how well it is possible to image and subtract 
out the point sources. 



We note that FIELDS I and IV have a relatively longer 
on-source observation time in comparison to FIELD II and 
HI, and this is reflected in the fact that FIELDS I and IV 
have a lower noise level and higher sensitivity in comparison 
to FIELD II and III (Table 2). The left panel of Figure 9 
shows a more detailed view of the brightest source for the 
most sensitive and the least sensitive fields among the four 
fields that we have imaged. For making these images, we 
have applied the appropriate phase shifts so as to bring the 
brightest source to the phase center of the image. The ratio 
of the peak flux density to the rms noise (Peak/Noise, Table 
2 ) has a maximum value of 700 in FIELD I, and a mini- 
mum value of (422) in FIELD III. The brightest source in 
each of our fields is also found to be accompanied by several 
regions of negative flux density. These are presumably the 
result of residual phase errors which were not corrected in 
our self calibration process. We find that the Peak negative 
flux density is 14mJy/Beam in FIELD I, and 47mJy/Beam 
in FIELD IV. These correspond to 1.5% and 5.3% of the 
Peak flux density in the respective FIELDS, the values lie 
within this range for the two other fields (Table 2). In ad- 
dition to the pixels with negative flux densities, we also see 
several regions of positive flux densities well above the 5a 
noise levels. Both the positive and negative regions are imag- 
ing artifacts which are possibly the outcome of calibration 
errors. The self-calibration steps implemented in the earlier 
stage of the analysis have considerably reduced the imag- 
ing artifacts. The artifacts, however, not entirely removed 
through self calibration. A visual inspection of the image of 
the brightest source in FIELD I (top left panel of Figure 9) 
shows that there are around 10 distinct features with flux 
densities around 7mJy arising from imaging artifacts. The 
artifacts effectively increase the local rms noise in the vicin- 
ity of the brightest source, and reduces the dynamic range 
of the image. 

We next consider how well we can subtract out the 
brightest source in the four fields that we have analyzed. 
We model the brightest source using the Clean Components 
of the continuum image. The visibilities corresponding to 
these clean Components were subtracted from the origi- 
nal full frequency resolution uv data using the AIPS task 
UVSUB. The right panel of Figure 9 shows the continuum 
image made with the residual visibilities after the brightest 
source was subtracted out. We find that the peak residual 
has values (—15, —64, —54, 34) mJy/Beam which correspond 
to (1.6,4.0,3.1,3.8)% of the original source that was sub- 
tracted out from the four fields respectively. FIELD I has 
the least residuals and the image is largely featureless in 
comparison to FIELD II which has the maximum residuals. 

The continuum image used for source subtraction does 
not take into account the possibility that the variation in the 
source flux across the observational band. In principle, we do 
not expect a significant flux variation across the relatively 
small frequency bands in our observation. Figure 10 shows 
the channel spectra through the brightest pixel in two of our 
observed fields. Contrary to our expectation, we find 15 % to 
30 % variations in the spectra. Further, these variations are 
not smooth and they show random and abrupt variations. 
The channels where data is missing have been flagged to 
avoid broad band RFI. This spectral variation is mainly due 
to errors in the bandpass calibration. It is very clear that 
it will not be very effective to model the observed spectral 



12 A. Ghosh, J. Prasad, S. Bharadwaj, S. S. Ali and J. N. Chengalur 



. 1 1 


i i 


1 1 T 


- - FIELD I 








■ ~ s 










- _<® rj 




V, - 




















\ J 






























o 






0, 


i i 


I 1 - 



05 29 45 



30 



15 00 28 45 

RIGHT ASCENSION (J2000) 
Cont peak flux = 8.9720E-01 JY BEAM 
Levs = 1.300E-03 • (-7, 7, 9, 15, 20, 50, 100, 300) 




05 29 45 30 15 00 28 45 

RIGHT ASCENSION (J2000) 
Cont peak flux = -1.581 8E-02 JY/BEAM 
Levs = 1.300E-03 * (-20, -15, -9, -7, 7, 9, 15, 
20) 




06 05 00 04 45 30 15 00 03 45 06 05 00 04 45 30 15 00 03 45 

RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 

Cont peak flux = 1.5670E+00 JY BEAM Cont peak flux = -6.4618E-02 JY BEAM 

Levs = 2.500E-03 ' (-20, -15, -9, -7, 7, 9, 15, Levs = 2.500E-03 * (-20, -15, -9, -7, 7, 9, 15, 

20,50,100,300,500,700) 20) 



Figure 9. This figure shows the position around the brightest source (marked with 'X') for FIELD I and II before (left panel) and after 
(right panel) the source is subtracted. 



variation of the source using a smooth polynomial or power 
law. It may be a better strategy to make separate images 
at each frequency channel and individually subtract out the 
CLEAN component from each channel of the uv data, pro- 
vided the signal to noise ratio is sufficient for single channel 
images. 

We next consider source subtraction from the entire 
field of view. Pixels with flux density above 7a times the rms 
noise were visually identified as sources and subtracted out 
using exactly the same procedure as used for the brightest 
source. It is expected that at this stage most of the genuine 
sources seen in Figure 1 have been removed from the uv data. 
Figure 11 shows the corresponding images made from the 
residual visibilities after source subtraction. The flux density 
in these images are in the range of (21, 39, 58, 30) mJy/Beam 
to (—14, — 29, —51, —51) mJy/Beam respectively for the four 



fields that we have analyzed. Most of the residuals, we find, 
are clustered in a few regions in the image. These regions 
have a typical angular extent of 15' and are centered on the 
locations of the bright sources in the FoV. The residuals are 
essentially imaging artifacts that were not modelled using 
Clean Components. The rest of the image, leaving aside a 
few isolated regions, is largely free of artifacts and devoid of 
any visible feature. We find that FIELD I has the highest 
sensitivity (lowest rms noise) amongst our observed fields. 
We also see that source removal is most effective for this 
field. We see (Figure 11) that most of the image is free of 
residual structures after source subtraction, except for two 
small regions where most of the artifacts are localised. It is 
quite evident that FIELD I is the best field and we use this 
for the entire subsequent analysis. 
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Figure 10. The figure shows the Flux/Beam vs channel (Frequency) plots at the brightest source position for two fields. Note that some 
channels were heavily corrupted by broad-band RFIs and we flagged those channels which are shown as 'gap' in FIELD I and II. 



5.1 Differential Source Count 

The differential source count plays a very important role in 
estimating the power spectrum of point sources. Further, 
this also provides invaluable information as to the nature of 
the radio sources. Here we evaluate the differential source 
count using the discrete sources identified in FIELD I down 
to a flux limit of 9 mjy. 

We use the AIPS task SAD to identify and extract the 
sources from the continuum image of FIELD I (Figure 1). 
SAD identifies potential sources based on the peak bright- 
ness and fits a Gaussian model to the sources. Source iden- 
tification was restricted within a radius of 1.45° from the 
phase centre, where the primary beam response is within 
50 per cent of its central value. We have used a conserva- 
tive peak brightness detection limit of 7a (> 9mJy/Beam 
) for the initial source selection using SAD. This minimizes 
the number of noise spikes that are spuriously detected as 
sources. 

As mentioned earlier, the local rms noise increases near 
the bright sources due to the presence of imaging artifacts. 
We have estimated the local rms in different parts of the 
residual image by applying the AIPS task RMSD. The local 
rms at each pixel was estimated using an area of approxi- 
mately 256 x 256 pixels centered on that pixel. Sources with 
a peak brightness in excess of 5a times the local rms noise 
were finally selected. Further, we have visually inspected the 
regions close to the brightest sources and discarded sources 
that appeared to be imaging artifacts. This ensures that we 
do not include any bright imaging artifacts as sources. The 
flux density of all the extracted sources were corrected for 
the GMRT primary beam using an eighth order polynomial 
(Kantharia & Pramesh Rao 2001). 

We have finally identified a total of 206 sources (Figure 
12) within a 1.45° radius from the pointing centre . The full 
source catalogue is given in the online version of the paper ( 
Appendix A) and is presented in order of increasing RA. For 
each source, the catalogue lists the position of the sources 
in RA, DEC, peak flux density, local rms noise at the source 
location, the integrated flux density and the error in the 
integrated flux density. 

We have used the source catalogue to determine the dif- 
ferential source count dN/dS, the number of sources dN in 
the flux interval dS. Note that S here refers to the integrated 



flux density of the source. To determine dN/dS we have 
binned the sources in the range 9.11 < S < 945.71 mjy into 
10 logarithmic bins in S (Table 3) and counted the number 
of sources (N) in each bin. We next consider two corrections 
that have to be considered in interpreting the source counts. 
The first is the Eddington bias where the random noise ar- 
tificially boosts the number of sources in the faintest bin. 
The 7a detection limit and the subsequent visual inspection 
help to avoid spurious sources and minimize the effect of 
the Eddington bias. Further, studies at 610 MHz, Moss et 
al. (2007) indicate that the Eddington bias only influences 
the faintest flux density bin, increasing the source count by 
approximately 20%. Since the number of sources in our low- 
est flux bin is comparatively small, we decided to make no 
correction for this effect. 

The second effect is related with the fact that extended 
sources with peak brightnesses below the survey limit but 
integrated flux densities above this limit would not be de- 
tected by our source detection procedure. It is possible to 
get an estimate of this effect with the knowledge of the an- 
gular size distribution of the sources as a function of flux 
density. This helps to identify the incompleteness due to ex- 
tended sources in the estimated source count distribution, 
well-known as resolution bias. Recently, Moss et al. (2007) 
have estimated that they will miss approximately 3% of the 
sources due to the resolution bias at 610 MHz. In our case, 
we do not find a significant number of extended sources and 
therefore we choose to make no correction for this effect. 

The actual threshold of flux density for source detection 
varies across the image, this being 5 times the local rms. 
Consequently, for each flux bin in Table 3 we only have a 
fraction / of the image where a source at the faint end of the 
bin can be detected. We correct for this by using N c — N / f, 
the noise corrected source count, to estimate the differential 
source count as 



dN 
dS 



AAS 



(8) 



where A is the total area of the image in steradians, 
and AS is the width of the flux bin in Jy. We have used 
the Poisson error AiV c = y/N c to estimate the error in the 
differential source count. 

We have used least square to fit a power law to the dif- 
ferential source count dN/dS as a function of the mean flux 
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Figure 11. This shows the same figure as the Figure 1 except that all the pixels with flux density above 7a noise value were visually 
identified as genuine sources and not artifacts have been fitted with clean components and removed from the visibility data from which 
this image was made. It is expected that most of the genuine sources have been removed from this data. 



density S of sources in each bin. We find the best fit power 
law to be dN/dS = l0 3 - 75±0 06 x S" 1 - 6 * 01 for our data. 
We note that Di Matteo et al. (2002) have found a single 
power-law fit dN/dS = 4000 x S" 1 ' 76 upto the flux den- 
sity of 880 mjy which is mostly consistent within the error 
bars of our fit. Figure 13 shows the normalised differential 
source count S 2 ' 5 x dN/dS as a function of the flux density 
S assuming an Euclidean Universe. We have compared our 
results with the 330 MHz source count model from Wieringa 
(1991) where the fit to the combined source count from six 
different fields (in the range from 4 mjy to 1 Jy) is given by 




= 0.976 + 0.61362 + 0.3028a; 2 - 0.083a: 3 (9) 



where x = log 10 (S) and S is given in mjy. 

This source count model was scaled down to our observ- 



ing frequency of 150 MHz assuming that the source fluxes 
can be scaled as S oc v~ a . We have tried out various val- 
ues of the mean spectral index a. We find that the model 
(eq. 9) is most consistent with our measurements for a mean 
spectral index of a — 0.7 (Figure 13). 

We find that our differential source count (Figure 13) 
which is determined over a flux range from 9 mjy to 1 Jy 
is well fitted by a single power-law of slope 0.97 ± 0.07. 
We note that George & Stevens (2008) have estimated a 
single power law of slope 0.72 in the flux range of 20mJy 
to 2 Jy from a 150 MHz GMRT survey centered around e 
Eridanus. In Figure 13 we show the normalised differential 
source count from their survey. We find that they have under 
determined the normalised differential source count in the 
entire flux range and it is most likely due to the low num- 
ber of sources (113) that were used to determine the source 
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Figure 12. The triangles show the 206 sources within a radius 
of 1.45° from the phase centre. Here, A x = COS(DECC)x(RA- 
RAC) and A y = (DEC-DECC) are the angular displacements 
from the phase center RAC and DECC. The RA and DEC rep- 
resents the positions of the different sources. 
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Figure 13. The figure shows the differential source counts at 
150 MHz (mc90 and ha88 (McGilchrist et al. 1990 and Hales, 
Baldwin & Warner 1988), George & Stevens (2008), Interna et al. 
(2011), Our Result), normalized by the value expected in a static 
Euclidean universe. A steady decrease in the source counts is 
observed as we approach lower flux values. The continuous curve 
represents the functional form of the fit (valid from 4 mjy to Uy ) 
to the source counts by Wieringa (1991), scaled down to 150 MHz 
assuming a mean spectral index of -0.70. 



Table 3. 150 MHz differential source counts for FIELD I. The 
columns show bin flux limits, the mean flux density of sources 
in each bin, number of sources, the noise corrected number of 
sources and dN/dS with error. 



Flux Bin 


5 


N 


N c 


dN c /dS 


(mjy) 


(mjy) 






(sr-ijy" 1 ) 


9.11 


- 14.49 


12.23 


9 


32.82 


(3.03 ± 0.53) x 10 6 


14.49 


- 23.06 


18.48 


25 


35.15 


(2.04 ± 0.34) x 10 6 


23.06 


- 36.68 


29.38 


11 


44.16 


(1.61 ± 0.24) x 10 6 


36.68 


- 58.35 


44.96 


38 


38.00 


(8.72 ± 1.41) x 10 5 


58.35 


- 92.82 


72.71 


27 


27.00 


(3.89 ± 0.75) x 10 5 


92.82 - 


- 147.67 


118.51 


23 


23.00 


(2.08 ± 0.43) x 10 5 


147.67 


- 234.91 


192.87 


18 


18.00 


(1.03 ±0.24) x 10 5 


234.91 


- 373.70 


306.16 


10 


10.00 


(3.58 ± 1.13) x 10 4 


373.70 


- 594.48 


495.95 


8 


8.00 


(1.80 ±0.64) x 10 4 


594.48 


- 945.71 


799.70 


1 


4.00 


(5.66 ±2.83) x 10 3 



count. For comparison we also plot the combined results 
from 6C and 7C survey 6 in Figure 13. We find that there 
is a good agreement between our source count and those 
derived from the 6C, 7C survey (Hales, Baldwin & Warner 
1988; McGilchrist et al. 1990) in the flux range 0.1 Jy to 
1.0 Jy. Interna et al. (2011) have determined the normalised 
differential source count at 153 MHz using a catalog of 598 
sources in the flux range of 4.75 mjy to 3 Jy. We notice that 
our source counts are roughly equal to their in the high flux 
range of 30.0 mjy to 1.0 Jy. Towards the lower flux ends 
(< 30.0 mjy) we find their source counts are increasingly 
higher compared to our source count. We note that our esti- 
mated source count agrees well with Ishwara- Chandra et al. 
(2010) in the low flux range (< 30.0 mjy) where they found 



http: / / web.oapd.inaf.it /rstools /sreent /150MHz.dat 



a relatively low population of sources compared to Interna 
et al. (2011). 

The simulated source count model proposed by Jack- 
son (2005) at 151 MHz predicts a flattening in the power 
law shape below 10 mjy due to a growing population of Fa- 
naroff - Riley (Fanaroff & Riley 1974) class I (FR I) radio 
galaxies. The turnover in the source count is a well-known 
observed feature at 1.4 GHz (Condon 1989; Hopkins et al. 
1998; Seymour et al. 2004) and occurs around 1 mjy, which 
is equivalent to ~ 5 mjy for an average spectral index of 0.7. 
We note that although this flattening is observed in deep 
radio surveys at higher frequencies (1.4 GHz, Windhorst et 
al. (1990); 610 MHz, Garn et al. 2007, 2008), our current 
survey depth is not sufficient to detect this flattening. This 
suggests that our catalog is dominated by the classical radio- 
loud AGN population which are the predominant sources at 
higher flux densities. 

We note that the predicted angular power spectrum Ce 
(Figure 6) does not change very significantly if we use the 
differential source count measured here instead of the fore- 
ground model used for the predictions in Section 4. 



6 DIFFUSE GALACTIC FOREGROUND 

The diffuse synchrotron radiation from our Galaxy domi- 
nates the sky at low radio frequencies like 150 MHz. Haslam 
et al. (1982) surveyed the full sky with an angular resolution 
of 0.85° at 408 MHz using single dish observations. The sky 
brightness temperature is found to vary in the range UK 
to 4, 247 K across the entire sky. Assuming that the syn- 
chrotron brightness temperature scales as T oc v~ 2 ' s with 
frequency (Platania et al. 2003), we have a conversion factor 
of 16.47 from 408 MHz to 150 MHz. This gives the brightness 
temperature range from 181 K to 69, 948 K across the en- 
tire sky at 150 MHz, and (660, 495, 330, 495) K at the respec- 
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tive phases centers of the four fields that we have observed 
(Table 1). The subsequent analysis is restricted to FIELD 
I where we have been able to achieve the best sensitivity, 
and point source subtraction is also most effective. Further 
this is the lowest Galactic latitude field (6 = 13.89°) among 
our observed fields and it may provide an upper limit on 
the expected diffuse foreground emission for four fields that 
we have observed. Figure 14 shows the 408 MHz (Haslam 
et al. 1982) brightness temperature distribution within a 
10° x 10° region of FIELD I. The figure also shows the 
150 MHz brightness temperature distribution obtained by 
scaling the 408 MHz maps. The scaling to 150 MHz takes 
into the angular variation of the spectral index (Platania et 
al. 2003), and we have used LFmap (E. Polisensky 2009) 
to generate the data for these figures. The smallest baseline 
in our observation is around U = 100, which corresponds 
to an angular scale of U^ 1 « 30'. Our observations are not 
sensitive to intensity variations at angular scales larger than 
this, and consequently we do not expect the structures seen 
in Figure 14 to be imprinted in our observation. Very little 
is known about the angular structure of the diffuse Galactic 
synchrotron radiation on sub-degree scales. 

Model predictions (Ali, Bharadwaj & Chengalur 2008), 
which extrapolate the statistical properties measured at 
large angular scales and higher frequencies, predict that 
point sources are the most dominant contribution at sub- 
degree scales. We expect the diffuse radiation to dominate 
if the point sources can be individually modelled and re- 
moved with high level of precession. Bernardi et al. (2009) 
have analyzed 150 MHz WSRT observations where they have 
subtracted out the point sources to reveal the structure of 
the fluctuations in the diffuse Galactic synchrotron radia- 
tion at 13' angular scales. We note that their observation 
was carried out at a low Galactic latitude (6 = 8°) where 
the Galactic emission is expected to be relatively larger than 
our targeted field. It is expected that we should be able to 
use the residual data to characterize the diffuse radiation 
provided the point sources have been removed to an ade- 
quate level of sensitivity. 

We have discussed point source subtraction in Section 5. 
For FIELD I, Figure 1 and Figure 11 show 4.0° x 4.0° contin- 
uum images of the field before and after source subtraction. 
All the pixels with flux density greater than lOmJy/Beam 
were visually inspected, and those which appeared to be 
genuine sources were fitted with Clean Components us- 
ing tight boxes. The continuum Clean Components were 
then subtracted from the visibility data. The image, after 
source subtraction, has residual flux density in the range of 
— 14 m Jy /Beam to 21 mJy/Beam arising from imaging arti- 
facts. We recollect here that these residuals are highly lo- 
calized in the vicinity of a few regions that contained the 
brightest sources (Figure 11). The bulk of the field is largely 
free of artifacts and is consistent with noise. 

It is possible that very bright sources beyond the 
4.0° x 4.0° field that we have imaged also contribute to the 
measured visibilities. To account for this possibility we have 
also imaged a 8.0° x 8.0° region after the sources within 
the central 4.0° x 4.0° region have been subtracted out. 
We have then removed all the sources automatically to a 
lOmJy/Beam level using the 8.0° x 8.0° image. 

We expect to see the diffuse background radiation after 
the point sources have been removed, however the high reso- 



lution (20" x 18") residual image does not exhibit any diffuse 
structure. The noise and the residuals after point source re- 
moval appear to dominate the high resolution image. We 
expect the Galactic synchrotron radiation to be relatively 
larger in comparison to the noise and residuals if we consider 
larger angular scales. We have considered the angular scale 
10', and following Bernardi et al. (2009) we have made two 
images which refer to this angular scale. Initially, we made 
an image which does not have any visibilities with baselines 
U < 170. This restriction imposes the condition that the 
resulting image does not contain any information on angular 
scales grater than 10'. For the second image (Figure 15), the 
visibilities were tapered with a Gaussian in the uv plane so 
as to produce a synthesised beam of FWHM 620" x 540". 
This image does not contain information at angular scales 
below 10'. 

We see that the image, made by including only the base- 
lines |U| > 170, which does not contain any information 
above 10' looks very similar to the high resolution image. 
Both of these images are dominated by the noise, and the 
residuals from the brightest point sources are seen to be lo- 
calized near the center of the image. In contrast, the low res- 
olution image (Figure 15), which does not contain informa- 
tion at angular scales below 10', shows fluctuations which are 
uncorrelated with the point source distribution seen in the 
high resolution image (Figure 1). The maximum and mini- 
mum values of flux density in Figure 15 are 113mJy/Beam 
and — 139mJy/Beam respectively, which are comparable to 
5cr where a (rms) is 23.5mJy/Beam (Table 4). The indi- 
vidual regions corresponding to the very high (or low) pixel 
values also have an angular extent that is comparable to the 
synthesized beam. We interpret these features, which cor- 
respond to brightness temperature fluctuations of the order 
of 20 K (Table 4) , as a tentative detection of the Galactic 
synchrotron radiation at the 10' angular scale. 

We next increase the uv taper to produce a synthesized 
beam with a larger FWHM of 1070" x 864" . Figure 16 shows 
the corresponding image which does not have any informa- 
tion at angular scales below approximately 16'. The maxi- 
mum and minimum values of flux density in this image are 
436 mjy /Beam and — 353mJy/Beam respectively, which are 
greater than 10a where a (rms) is 35mJy/Beam (Table 4). 
The individual regions corresponding to the very high (or 
low) pixel values have an angular extent that is compara- 
ble to, if not bigger than, the synthesized beam. In fact the 
size of these regions approaches the largest angular scales 
(~ 30') that can be probed in our observation. We interpret 
these lOcr fluctuations as a detection of the fluctuations in 
the Galactic synchrotron radiation at the 16' angular scale. 
Converting to brightness temperature (Table 4), we have a 
peak fluctuation of 26 K in our field of view. 

Based on our observation we conclude that the residu- 
als after point source subtraction represent structures in the 
Galactic synchrotron radiation on angular scales of 10' to 
20'. The noise and the residuals from point source subtrac- 
tion, however, dominate the data at smaller angular scales 
where we are unable to make out the Galactic synchrotron 
radiation. 
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Figure 14. The figure in the left panel shows the brightness temperature distribution around FIELD I at 408 MHz (Haslam ct al. 1982) 
after destriping (Platania et al. 2003). In the right panel we show the temperature values at 150 MHz where we scale the temperature 
values with the average spectral index map from (Platania et al. 2003). We use (E. Poliscnsky 2009) to generate this figures. 
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Figure 15. This shows the image where we tapered the uv plane at |U| = 200. The synthesized beam has a FWHM of 620" X 540". 
All the genuine point sources were removed down to lOmJy level. Diffuse structures begin to appear on > 10 scales. The brightest 
structures in this map are at a ~ 5u level compared to the local rms of 23.5 mJy/Bcam. 
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Figure 16. This shows the image where we tapered the uv plane at |U| = 100. The synthesized beam has a FWHM of 1070" X 864". 
All the genuine point sources were removed down to lOmJy level. The brightest structures in this map are at a 10ct level compared to 
the local rms of 35mJy/Beam. 



Table 4. Rms fluctuations in the FIELD I as a function of angular 
resolution 



Angular 


rms 


Conversion factor 


rms 


resolution 


(mJy/Beam) 


(mJy/Beam) to (K) 


(K) 


20" x 18" 


1.30 


171.80 


223.34 


620" x 540" 


23.50 


0.17 


4.00 


1070" x 864" 


35.00 


0.06 


2.1 



6.1 Power Spectrum 

We have used the visibilities, after point source subtraction, 
to estimate the angular power spectrum Ce = Ce(Au = 0). 
The angular power spectrum Ce (Figure 17) clearly shows 
two different scaling behavior as a function of I. At low angu- 
lar multipoles (£ < 800), which correspond to angular scales 
larger than 10', we find a steep power law behavior which 



is typical of the Galactic synchrotron emission observed at 
higher frequencies and larger angular scales (Bennett et al. 
2003). The angular power spectrum flattens out for I > 800, 
and we find that it remains nearly flat out to £ ^ 8, 000 
which corresponds to angular scales of ~ 1 . The nearly flat 
region arises from the point sources which have a flux that 
is below the threshold for source identification and removal, 
and hence have not been subtracted from the data. The er- 
ror's in modelling and subtracting the identified sources also 
contribute to this. 

It is clear from this analysis that Ce, after point 
source subtraction, is dominated by the diffuse Galactic syn- 
chrotron emission at I < 800 which corresponds to ~ 10'. 
We also note that the angular power spectrum Ce was di- 
rectly estimated from the visibilities using the method dis- 
cussed in Section 3 which makes absolutely no reference to 
the image. However, it is reassuring to note that the mea- 
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sured Ce is quite consistent with the behavior that we no- 
ticed in the image plane. We do not find any significant 
features, other than the imaging artifacts, in the high reso- 
lution residual image and in the residual image made by dis- 
carding the short baselines (|Uj > 170). Significant features 
only appear in the residual image if the angular resolution 
is > 10'. The measured angular power spectrum predicts 
the f luctuations in the Galactic synchrotron emission to be 
~ y / £(£ + l)C e /2iv ~ 10 K at £ = 800 which is comparable 
to the 20 K features seen in the 10' resolution image (Figure 
15). We have used a weighted least square to fit a power law 

Ct = Ax(^y m K> (10) 

to the measured Ce for £ < 800 . We find the best fit 
values A = 513 ± 41 mK 2 and /3 = 2.34 ± 0.28 for which Ce 
is shown in Figure 17. 

Bernardi et al. (2009) have analyzed 150 MHz WSRT 
observations, where they have subtracted out point sources 
above 15 mjy and used the resulting image to estimate the 
angular power spectrum Ce- We note that their analysis dif- 
fers from our in that they have estimated Ce from the image 
whereas we have directly used the measured visibilities to es- 
timate Ce (Figure 17). Our findings, however, are very simi- 
lar to those of Bernardi ct al. (2009) who find that the mea- 
sured Ce is well described by a power law for £ < 900. Their 
best fit parameters are A = 253±40mK 2 and /3 = 2.2 ±0.3. 
The slope, we note, is consistent with our findings, the am- 
plitude, however, is half our value. The difference in ampli- 
tude is not surprising as the two different values refer to 
observations in two different parts of the sky. The ampli- 
tude measured by Bernardi et al. (2009) refers to a field 
with Galactic latitude b = 8°, whereas FIELD I of our ob- 
servations is at b = 13.89°. We would generally expect lesser 
synchrotron emission in our field which is at a higher ele- 
vation from the Galactic plane, however, the amplitude of 
Ce shows the opposite behaviour. The variation of the sub- 
degree scale angular structure of the Galactic synchrotron 
emission across different parts of the sky is not known at 
present, and it is not possible to make any conclusive state- 
ment about this apparent discrepancy. It is possible that 
there may be some extra power coming from the clustering 
of unresolved point sources below I < 800. We note that 
our estimate of the angular power spectrum for the Galac- 
tic diffuse emission (Figure 17) matches quite well with the 
foreground model prediction of Ali, Bharadwaj & Chengalur 
(2008), assuming that all point sources above a threshold 
flux limit of 20 mjy (see for reference top left panel of Fig- 
ure 11, FIELD I) have been removed. 

Strictly speaking, it is not justified to compare the am- 
plitude of the angular power spectrum obtained in observa- 
tions at two different frequencies unless their slopes (an- 
gular spectral index) are consistent with each other. We 
have not found any other observations (except Bernardi et 
al. (2009)) of the northern Galactic plane where the angu- 
lar spectral index is consistent with our findings. Despite 
this, we have used an earlier work (La Porta et al. 2008) 
at higher frequencies to estimate the amplitude of the an- 
gular power spectrum expected at our observing frequency. 
La Porta et al. (2008) have estimated angular power spec- 
trum of all-sky total intensity maps at 408 MHz (Haslam 
et al. 1982) and 1420 MHz (Reich 1982; Reich & Reich 



1986; Reich et al. 2001). The angular power spectrum of 
the Galactic synchrotron emission is measured down to the 
angular mode of £ = 200 and I = 300 at 408 MHz and 1420 
MHz respectively. Using the best fit parameters (tabulated 
at £ = 100) at 408 MHz and 1420 MHz, we obtain the am- 
plitude of the angular power spectrum Ce=200 — 230 mK 2 
and 0.22 mK 2 respectively close to the Galactic latitude of 
b = 15° (close to FIELD I). By extrapolating the values 
of Ce=2oo from 408 MHz and 1420 MHz using a mean fre- 
quency spectral index of a = 2.5 (de Oliveira-Costa et al. 
2008) (Ce oc v~ 2a ) we find that the expected contribu- 
tions at 150 MHz from the Galactic synchrotron emission 
are 3.42 x 10 4 mK 2 and 2.08 x 10 4 mK 2 respectively. In our 
observation we find Ce=2m = 2.22 x 10 4 mK 2 which lies 
within the range 2.08 x 10 4 mK 2 to 3.42 x 10 4 mK 2 , and 
hence we conclude that our finding is consistent with the 
values extrapolated from higher frequency and larger angu- 
lar scales. 

Giardino et al. (2001) have analyzed the fluctuations 
in the Galactic synchrotron radiation using the 2.3 GHz 
Rhodes Survey where they find a slope /3 — 2.43 ± 0.01 
across the range 2 < I < 100. Giardino et al. (2002) have 
analyzed the 2.4 GHz Parkes radio continuum and polariza- 
tion survey of the southern Galactic plane where they find 
a slope p = 2.37 ± 0.21 across the range 40 < I < 250. 
Our slope, measured at smaller angular scales, is consistent 
with these findings. There has been considerable work on 
modelling the Galactic synchrotron radiation at the higher 
frequencies relevant for the Cosmic microwave background 
radiation (CMBR) (tens of GHz). These models predict a 
power law scaling behaviour Ce oc £~' 9 where /3 have values 
in the range 2.4 to 3 down to £ = 900 (Tegmark & Efstathiou 
1996; Tegmark et al. 2000). 

We note that the slope measured by us at 150 MHz 
is consistent with these model predictions. (Bennett et al. 
2003) have determined the angular power spectrum of the 
Galactic synchrotron radiation using the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) data in the frequency 
range 23 to 94 GHz. Their measurement is restricted to 
£ < 200 where they find a scaling Ce ~ £~ 2 which is slightly 
flatter than the slope we obtain at smaller angular scales. 



7 DISCUSSION AND CONCLUSIONS 

We have analyzed 150 MHz GMRT observations in four dif- 
ferent fields of view (Table 1, Figure 1). We have used the 
multi- frequency angular power spectrum (MAPS) Ce(Av) 
to jointly characterize the statistical properties of the fluc- 
tuations as a function of the angular multipole £ and the 
frequency separation Av (Figure 4) across the range 700 < 
£ < 2 x 10 4 and < Av < 2.5 MHz. We find that the mea- 
sured Ce has values around 10 4 mK 2 at £ ~ 1000, and it 
drops by around 50% at £ ~ 10 4 . The measured Ce is fore- 
ground dominated and is more than 7 orders of magnitude 
larger than the expected HI signal. 

Analytic estimates of the HI signal (Bharadwaj and 
Ali 2005) show that we expect C e HI (Av) to decorrelate 
rapidly with increasing Av, and fall by 90% or more at 
Av > 0.5 MHz . In contrast, the foregrounds from differ- 
ent astrophysical sources are expected to be correlated over 
large frequency separations (Av > 2 MHz). This holds the 
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Figure 17. The measured Ce before (triangle), and after (circle) 
source subtraction, with l<r error-bars shown only for the latter. 
The dashed line shows the best fit power-law (I < 800) after 
source removal. The thin solid line (lowest curve) with 1 - a error- 
bars shows the model prediction of Ali, Bharadwaj & Chengalur 
2008 assuming that all point sources above S c = 20 mjy have 
been removed. 

promise of allowing us to separate the signal from the fore- 
grounds (Ghosh et al. 2011b). We find that the measured 
Ce(Av) (for reference see Figures 7 and 8) shows a smooth 
Av dependence with a 20 to 40 % variation across the mea- 
sured 2.5 MHz Av range. However, contrary to our expec- 
tations, in addition to a smooth behaviour we also notice 
abrupt variations and oscillations (< 10%) in the Av de- 
pendence of the measured Ct(Av). These abrupt variations 
and oscillations, whose origin is at present not understood, 
poses a severe impediment for foreground removal. Polyno- 
mial fitting along the frequency axis has been extensively 
considered for foreground removal. All attempts in this di- 
rection may fail due to the presence of abrupt variations 
and small oscillations which cannot be filtered out with low 
order polynomials. 

The measured Ce's, we find, are consistent with fore- 
ground models which predict that extragalactic point 
sources make the most dominant contribution at the sub- 
degree scales that we have probed here (Figure 6). The 
brightest source (~ 1 Jy) in each field alone contributes 
around 10% of the total measured Ce- It is very important 
to correctly identify the point sources and subtract these 
out at a high level of precision. We have carried out point 
source subtraction in all the four fields that we have ana- 
lyzed here (Figures 9 and 11). Considering only the bright- 
est source, we find that source subtraction is most effective 
in FIELD I, where the image also has the lowest rms noise 
of 1.3mJy/Beam (Table 2). The residual artifacts after sub- 
tracting out the brightest source, we find, are within 1.5% 
of the of the original source flux (905 mjy). 

The subsequent analysis was restricted to FIELD I 
which has the lowest rms noise. The Clean Components cor- 
responding to all the sources above 10 mjy were subtracted 
out from the visibility data using the AIPS task UVSUB. We 
find that the resulting angular power spectrum Ce falls to 
~ 1, 000 mK 2 in the I range 800 to 8, 000 (Figure 17), which 
is roughly one-tenth of the Ce before source subtraction. At 
these multipoles, we interpret the measured Ce, after source 



subtraction, as arising from a combination of the resid- 
ual artifacts from the bright sources and the faint sources 
(< lOmJy) that have not been removed. The behaviour at 
lower multipoles (£ < 800), we find, is well fitted by a power 
law C e = (513 ± 41) x (1000/£) 234±0 - 28 mK 2 which we in- 
terpret as the contribution from the diffuse Galactic syn- 
chrotron radiation. The measured slope is consistent with 
earlier WSRT 150 MHz observations (Bernardi et al. 2009), 
and also with 2.3 and 2.4 GHz results at smaller £ (Gia- 
rdino et al. 2001, 2002), whereas WMAP finds a flatter slope 
Ce ~ £~ 2 (Bennett et al. 2003) at smaller £ and much higher 
frequencies (23 - 94 GHz). 



8 ACKNOWLEDGMENT 

AG would like to thank M.H.Wieringa for sharing the 3rd 
chapter of his thesis. AG thanks Emil Polisensky for a de- 
tail discussion on LFmap. AG acknowledges CSIR, India for 
providing financial support through Senior Research Fellow- 
ship (SRF). SSA acknowledges the support by DST, India, 
under project No.: SR/FTP/PS-088/2010. The data used in 
this paper were obtained using GMRT. The GMRT is run 
by the National Centre for Radio Astrophysics of the Tata 
Institute of Fundamental Research. We thank the GMRT 
staff for making these observations possible. 



REFERENCES 

Ali S. S., Bharadwaj S.,& Chengalur J. N., 2008, MNRAS, 
385, 2166A 

Bharadwaj S. , & Ali S. S. 2005, MNRAS, 356, 1519 
Becker,R.H.,et al.,2001,AJ, 122, 2850 

Begum, A., Chengalur, J. N., & Bhardwaj, S. 2006, MN- 
RAS, 372, L33 

Bennett, A. S. 1962, MNRAS, 125, 75 

Bennett C.L., Hill R.S., Hinshaw. G. et al., 2003, ApJS, 
148, 97 

Bernardi, G., et al. 2009, A&A, 500, 965 
Bernardi, G., Mitchell, D. A., Ord, S. M., et al. 2011, MN- 
RAS, 413, 411 

Bowman J.D., Morales M.F. & Hewitt J.N., 2009, ApJ, 
695, 183 

Bharadwaj, S. & Sethi,S. 2001, JApA, 22, 293 
Chapman, E., Abdalla, F. B., Harker, G., et al. 2012, MN- 
RAS, 423, 2518 
Condon, J. J. 1989, Ap.J, 338, 13 

Cooray, A., & Furlanetto, S. R. 2005, MNRAS, 359, L47 
de Oliveira-Costa, A., Tegmark, M., Gaensler, B. M., et al. 

2008, MNRAS, 388, 247 

Di Matteo, T., Perna, R., Abel, T. & Rees, M.J., 2002, 

Ap.J, 564, 576 
Dunlop, J. S., & Peacock, J. A. 1990, MNRAS, 247, 19 
Dutta P., Begum A., Bharadwaj S., &: Chengalur J. N. 

2009, MNRAS, 398, 887 

Edge, D. O., Shakeshaft, J. R., McAdam, W. B., Baldwin, 
J. E., & Archer, S. 1959, Mem. R. Astron. Soc, 68, 37 

E. Polisensky, 'LFmap: A Low Frequency Sky 
Map Generating Program', Long Wavelength Ar- 
ray Memo Series No. Ill, Sep 7, 2007. [online] 
http://www.phys.unm.edu/lwa/memos. 



Foreground for EoR : 150 MHz GMRT observation 21 



Fan, X., ct al. 2002, AJ, 123, 1247 
Fanaroff, B. L., & Riley, J. M. 1974, MNRAS, 167, 31P 
Furlanctto , S. R. , Oh ,S. P.,. & Briggs,F., 2006, Phys.Rept. 
433, 181 

Gleser, L., Nusser, A., & Benson, A. J. 2008, MNRAS, 391, 
383 

Garn T., Green D.A., Hales S.E.G., Riley J.M. & Alexan- 
der P., 2007, MNRAS, 376, 1251 

Garn, T., Green, D. A., Riley, J. M., & Alexander, P. 2008, 
MNRAS, 387, 1037 

Ghosh, A., Bharadwaj, S., Ali, S. S., & Chengalur, J. N. 
2011a, MNRAS, 411, 2426 

Ghosh, A., Bharadwaj, S., Saiyad Ali, S., & Chengalur, 
J. N. 2011b,MNRAS, 418, 2584 

Giardino, G., Banday, A. J., Fosalba, P., et al. 2001, A&A, 
371, 708 

Giardino, G., Banday, A. J., Gorski, K. M., et al. 2002, 
A&A, 387, 82 

Ginzburg, V. L. & Syrovatskii, S. I., 1969, Ann.Rev.Astron. 

Astrophys., 7, 375 
George, S. J., & Stevens, I. R. 2008, MNRAS, 390, 741 
Hales, S. E. G., Baldwin, J. E., & Warner, P. J. 1988, MN- 
RAS, 234, 919 

Hales, S. E. G., Riley, J. M., Waldram, E. M., Warner, 
P. J., & Baldwin, J. E. 2007, MNRAS, 382, 1639 

Harker, G., et al. 2009, MNRAS, 397, 1138 

Harker, G., Zaroubi, S., Bernardi, G., et al. 2010, MNRAS, 
405, 2492 

Haslam, C. G. T., Salter, C. J, Stoffel, H. & Wilson, W. 
E., 1982, A&AS, 47, 1. 

Hopkins A. M., Mobasher B., Cram L., & Rowan-Robinson 
M. 1998, MNRAS, 296, 839 

Interna, H. T., van Wceren, R. J., Rottgering, H. J. A., & 
Lai, D. V. 2011, A&A, 535, A38 

Ishwara-Chandra, C. H., Sirothia, S. K., Wadadekar, Y., 
Pal, S., & Windhorst, R. 2010, MNRAS, 405, 436 

Ishwara-Chandra, C. H., Sirothia, S. K., Wadadekar, Y., 
Pal, S., & Windhorst, R. 2011, VizieR Online Data Cat- 
alog, 740, 50436 

Jackson, C. A., & Wall, J. V. 1999, MNRAS, 304, 160 

Jackson, C. 2005, PASA, 22, 36 

Jonas, J. L., Baart, E. E., & Nicolson, G. D. 1998, MNRAS, 
297, 977 

Jelic, V., Zaroubi, S., Labropoulos, P., et al. 2008, MNRAS, 
389, 1319 

Datta K. K., Roy Choudhury, T.,& Bharadwaj. S, 2007, 

MNRAS, 378, 119 
Kantharia N., Pramesh Rao A., 2001, GMRT Technical 

Note R00185 

Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, Ap.JS, 
180, 330 

Laing, R. A., Riley, J. M. & Longair, M. S. 1983, MNRAS, 
204, 151 

La Porta, L., Burigana, C, Reich, W., & Reich, P. 2008, 

A&A, 479, 641 
Liu A., Tegmark M., Bowman J.D., Hewitt J.N. & Zaldar- 

riaga M., 2009b, MNRAS, 398, 401 
Liu A., Tegmark M., & Zaldarriaga M. 2009, MNRAS, 394, 

1575 

McGilchrist, M. M., Baldwin, J. E., Riley, J. M., Titter- 
ington, D. J., Waldram, E. M., & Warner, P. J. 1990, 
MNRAS, 246, 110 



Morales, M. F., & Wyithe, J. S. B. 2010, ARA&A, 48, 127 
Moss D., Seymour N., McHardy I.M., Dwelly T., Page M.J. 

& Loaring N.S., 2007, MNRAS, 378, 995 
Oort M.J.A., 1988, A&A, 193, 5 

Peacock, J. A. 1999, Cosmological Physics, by John A. Pea- 
cock, pp. 704. ISBN 052141072X. Cambridge, UK: Cam- 
bridge University Press, January 1999 

Pen, U.-L., Chang, T.-C, Hirata, C. M., et al. 2009, MN- 
RAS, 399, 181 

Perley, R.A. 1999, ASP Conference Series, "Synthesis 
Imaging in Radio Astronomy II", Eds. G. B. Taylor, C. 
L. Carilli, and R. A. Perley, Vol. 180, p. 19 
Dutta, P. 2011, arXiv:1102.4419 
Petrovic, N., & Oh, S. P. 2011, MNRAS, 292 
Pindor, B., Wyithe, J. S. B., Mitchell, D. A., et al. 2011, 
PASA, 28, 46 

Platania, P., Burigana, C, Maino, D., et al. 2003, A&A, 
410, 847 

Platania, P., Bensadoun, M., Bersanelli, M., de Amici, G., 
Kogut, A., Levin, S., Maino, D., & Smoot, G. F. 1998, 
Ap.J, 505, 473 

Prasad, J., & Chengalur, J. 2012, Experimental Astronomy, 
33, 157 

Reich, W., 1982, A&AS, 48, 219. 
Reich, P., & Reich, W. 1986, A&AS, 63, 205 
Reich, P., Testori, J. C, & Reich, W. 2001, A&A, 376, 861 
Reich, P. & Reich, W., 1988, A&AS, 74, 7. 
Santos M. G., Cooray A., & Knox L. 2005, Ap.J, 625, 575 
Seymour, N., McHardy, I. M., & Gunn, K. F. 2004, MN- 
RAS, 352, 131 

Shaver, P. A., Windhorst, R, A., Madau, P.& de Bruyn, A. 
G.,1999,Astron. & Astrophys. ,345,380 

Sirothia, S. K., Saikia, D. J., Ishwara-Chandra, C. H., & 
Kantharia, N. G. 2009, MNRAS, 392, 1403 

Sun, X. H., Reich, W., Waelkens, A., & Enfilin, T. A. 2008, 
A&A, 477, 573 

Swarup, G., Ananthakrishnan, S., Kapahi, V. K., Rao, 
A. P., Subrahmanya, C. R., & Kulkarni, V. K. 1991, CUR- 
RENT SCIENCE, 60, 95 

Tegmark M. & Efstathiou G., 1996, MNRAS, 281, 1297 

Tegmark, M., Eisenstein, D. J., Hu, W., & de Oliveira- 
Costa, A. 2000, Ap.J, 530, 133 

Waldram, E. M. 1998, Observational Cosmology with the 
New Radio Surveys, 226, 63 

Windhorst R., Doug M. & Lyman N., 1990, ASPC, 10, 389 

Wieringa, M.H., 1991, PhD Thesis, Leiden University, Lei- 
den, The Netherlands 

Zaroubi, S., de Bruyn, A. G., Harker, G., et al. 2012, 
arXiv: 1205.3449 

Zaldarriaga, M., Furlanetto, S. R., & Hernquist, L. 2004, 
Ap.J, 608, 622 



APPENDIX A: SOURCE CATALOGUE 
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Table Al. The complete 150 MHz GMRT source catalogue, listed in ascending order of right ascension. The full table is available online. 



Source Name 


RA 
J2000.0 


Dec 
J2000.0 


Peak 
mjy beam -1 


Local Noise 
mjy beam -1 


Int. Flux Density 
mjy 


Error 
mjy 


GMRTJ051850.9+601311.6 


05: 


18:50.91 


+60:13:11.57 


43.74 


2.28 


66.22 


5.92 


GMRTJ051931. 0+601053.1 


05: 


:19:30.96 


+60:10:53.06 


37.39 


2.04 


38.08 


4.22 


GMRTJ051933.8+603114.4 


05: 


19:33.75 


+60:31:14.40 


31.42 


2.38 


37.90 


4.31 


GMRTJ051946.6+603205.3 


05' 


:19:46.62 


+60:32:05.28 


125.07 


2.34 


206.97 


6.11 


GMRTJ051958.5+600249.9 


05 


19:58.52 


+60:02:49.91 


18.84 


1.91 


19.08 


3.68 


GMRTJ052010.1+603703.0 


05 


20:10.06 


+60:37:03.02 


25.16 


2.22 


29.63 


4.86 


GMRTJ052016.0+601339.4 


05: 


:20:16.04 


+60:13:39.45 


31.43 


1.92 


42.68 


4.61 


GMRTJ052052.9+600426.4 


05: 


:20:52.94 


+60:04:26.38 


108.84 


1.86 


128.28 


3.87 


GMRTJ052059.2+593306.4 


05: 


:20:59.24 


+59:33:06.42 


23.36 


2.02 


37.92 


5.07 


GMRTJ052146.0+605705.9 


05 


:21:46.01 


+60:57:05.94 


47.74 


2.14 


68.69 


5.69 


GMRTJ052147.6+591732.7 


05: 


21:47.63 


+59:17:32.66 


152.35 


2.06 


226.12 


4.94 


GMRTJ052150.1+603309.9 


05: 


21:50.06 


+60:33:09.86 


21.56 


1.90 


22.62 


3.67 


GMRTJ052203.1+603652.6 


05: 


:22:03.13 


+60:36:52.65 


30.52 


1.90 


34.02 


3.21 


GMRTJ052207.7+594823.6 


05: 


:22:07.70 


+59:48:23.63 


18.81 


2.01 


24.49 


3.65 


GMRTJ052207.7+593820.6 


05 


:22:07.71 


+59:38:20.57 


35.74 


1.92 


39.56 


3.17 


GMRTJ052211.4+590248.1 


05: 


:22:11.40 


+59:02:48.13 


30.38 


2.39 


33.41 


4.64 


GMRTJ052212.5+602833.9 


05: 


:22:12.54 


+60:28:33.94 


17.41 


1.85 


30.61 


4.75 


GMRTJ052213.2+595938.8 


05 


:22:13.21 


+59:59:38.81 


156.82 


1.82 


166.62 


3.16 


GMRTJ052214.0+593104.8 


05 


22:13.98 


+59:31:04.78 


19.71 


1.88 


48.47 


6.01 


GMRTJ052214.9+590800.9 


05 


:22:14.89 


+59:08:00.94 


25.41 


2.24 


28.28 


4.37 


GMRTJ052237.8+594709.3 


05: 


:22:37.79 


+59:47:09.29 


14.87 


2.00 


14.99 


3.00 


GMRTJ052244.9+595140.3 


05: 


:22:44.95 


+59:51:40.31 


17.07 


1.93 


37.86 


4.86 


GMRTJ052246.5+591939.2 


05: 


:22:46.50 


+59:19:39.15 


43.15 


1.92 


49.83 


3.78 


GMRTJ052255.4+595129.7 


05 


:22:55.44 


+59:51:29.70 


472.06 


1.93 


684.06 


3.60 


GMRTJ052308.4+603902.6 


05' 


:23:08.44 


+60:39:02.55 


73.60 


1.74 


84.46 


3.60 


GMRT J052311. 8+600747.8 


05: 


23:11.75 


+60:07:47.84 


98.09 


1.73 


101.68 


2.87 


GMRTJ052317.8+595149.3 


05 


:23:17.84 


+59:51:49.26 


25.85 


1.92 


51.16 


4.27 


GMRTJ052320.8+601304.0 


05. 


:23:20.77 


+60:13:03.98 


492.90 


1.84 


552.31 


3.00 


GMRTJ052321.4+594929.7 


05: 


:23:21.42 


+59:49:29.69 


19.35 


1.97 


32.68 


3.84 


GMRTJ052324.6+594703.6 


05' 


:23:24.62 


+59:47:03.62 


275.72 


1.97 


566.23 


4.39 


GMRTJ052326.4+592542.5 


05 


:23:26.41 


+59:25:42.55 


15.58 


1.75 


17.61 


3.00 


GMRTJ052326.8+590426.4 


05 


:23:26.80 


+59:04:26.43 


109.28 


2.05 


110.68 


3.88 


GMRTJ052337.4+591059.7 


05: 


:23:37.40 


+59:10:59.68 


53.94 


1.95 


59.38 


3.35 


GMRTJ052337.8+600018.7 


05: 


23:37.83 


+60:00:18.68 


23.14 


1.66 


24.24 


2.62 


GMRTJ052356.3+603833.7 


05 


:23:56.29 


+60:38:33.69 


60.73 


1.65 


110.48 


4.52 


GMRTJ052402.5+593800.5 


05 


:24:02.55 


+59:38:00.52 


15.94 


1.63 


46.61 


5.60 


GMRTJ052408.9+585659.4 


05: 


:24:08.87 


+58:56:59.40 


205.90 


2.12 


225.58 


4.24 


GMRTJ052412.9+594652.1 


05: 


24:12.91 


+59:46:52.06 


296.69 


1.77 


338.91 


2.83 


GMRTJ052414.2+600001.3 


05 


: 24:14.21 


+60:00:01.26 


79.39 


1.62 


83.19 


2.64 


GMRTJ052416.1+601830.4 


05: 


24:16.10 


+60:18:30.44 


14.86 


1.75 


15.66 


2.75 


GMRTJ052417.4+602737.2 


05: 


:24:17.39 


+60:27:37.23 


46.03 


1.57 


52.44 


3.01 


GMRTJ052418.4+602708.0 


05' 


24:18.37 


+60:27:08.00 


55.73 


1.57 


57.65 


2.68 


GMRTJ052424.9+592200.1 


05 


:24:24.90 


+59:22:00.14 


16.30 


1.70 


18.10 


2.76 


GMRTJ052426.9+595234.0 


05: 


:24:26.94 


+59:52:34.00 


21.92 


1.69 


24.71 


2.73 


GMRTJ052433.3+591759.1 


05: 


:24:33.28 


+59:17:59.13 


23.50 


1.75 


34.46 


3.86 


GMRTJ052436.8+610912.2 


05: 


:24:36.78 


+61:09:12.23 


32.26 


1.99 


36.01 


4.55 


GMRTJ052440.6+585854.4 


05 


:24:40.56 


+58:58:54.44 


126.28 


2.07 


146.51 


4.12 


GJV1R1 J052448. 1+602715.1 


05: 


: 24:48.11 


+60:27:15.11 


13.49 


1.56 


14.77 


2.60 


GMRTJ052452.7+595004.5 


05' 


:24:52.67 


+59:50:04.49 


16.56 


1.67 


21.46 


2.90 


GMRTJ052452.9+605641.9 


05: 


:24:52.91 


+60:56:41.90 


30.12 


1.78 


31.98 


3.26 


GMRTJ052456.7+584255.2 


05 


:24:56.75 


+58:42:55.18 


42.00 


2.42 


43.78 


4.84 


GMRTJ052506.4+600128.6 


05 


:25:06.39 


+60:01:28.62 


14.41 


1.53 


15.20 


2.29 


GMRTJ052529.0+605027.3 


05' 


:25:29.01 


+60:50:27.33 


56.66 


1.62 


75.25 


3.75 


GMRTJ052533.3+603155.1 


05: 


:25:33.27 


+60:31:55.06 


134.08 


1.55 


176.86 


3.13 


GMRTJ052536.5+594516.9 


05: 


:25:36.48 


+59:45:16.90 


11.05 


1.68 


12.30 


2.31 


GMRT J052541. 5+601245.2 


05: 


:25:41.46 


+60:12:45.20 


14.55 


1.70 


34.63 


4.16 


GMRTJ052548.4+591802.6 


05: 


:25:48.41 


+59:18:02.63 


67.43 


1.66 


72.51 


2.96 


GMRTJ052600.1+603314.9 


05: 


26:00.10 


+60:33:14.93 


13.35 


1.52 


16.32 


2.94 


GMRTJ052645.5+583921.0 


05: 


:26:45.51 


+58:39:20.99 


80.73 


2.36 


83.13 


4.53 


GMRTJ052647.5+601450.0 


05 


:26:47.48 


+60:14:49.97 


47.41 


1.72 


200.70 


6.16 



